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The noise introduced into charge packets transferred through and 
stored in charge-transfer devices is calculated in a manner that includes 
all important relaxation, suppression, and correlation effects. First, the 
noise induced into each packet during each transfer phase from thermal, 
trapping, emission-current, and leakage-current fluctuations, whose 
statistics are nonstationary , and from clock-voltage fluctuations, whose 
statistics are stationary, is determined. Relaxation of the transferring 
charge to these fluctuations is found to suppress their size. Second, the 
accumulation (collecting) of the noise as each packet is transferred through 
the device is calculated neglecting the role of incomplete charge transfer. 
Attention is drawn to the significant differences between the collecting of 
storage-process noise, which is unsuppressed, transfer-process noise, 
whose spectral density is nearly totally suppressed at low frequencies, and 
modulation noise, which is nearly totally suppressed for digital and 
analog signals. Third, the role of incomplete charge transfer in suppressing 
the collecting of the noise is shown for digital signals and indicated for 
analog signals. We conclude with a numerical calculation of the maximum 
possible signal-to-noise ratio that can be expected from charge-transfer de- 
vices. The presentation is sufficiently general and detailed that, with a 
minimum of background informal noise theory, one can use the approach 
to evaluate noise in many novel, solid-state devices. 

I. INTRODUCTION 

Charge-transfer devices (CTD's), such as the bucket-brigade 1 (BB) 
and charge-coupled-device 2 (CCD) shift registers, are currently of 
great interest. These devices consist of a chain of charge-storage 
elements along which charge packets are transferred from input to 
output. Noise accompanying each individual transfer of each charge 
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packet will be introduced into each packet. It is the purpose of this 
article to calculate this noise and its cumulative effects on the output 
signal. 

Noise generated in solid-state devices has been treated extensively, 8-6 
and this prior work will be of great assistance to us here. However, 
there are three major, significant differences between much of this 
prior work and the present treatment. 7 ' 8 The first difference arises 
from the nature of the transfer process. 9 The usual treatments of 
noise discuss situations in which the (noiseless) currents, charge 
densities, conductances, etc. associated with the signal are time- 
invariant. Under such conditions, the statistics of the noise are also 
time-invariant and the noise is said to be stationary. Stationary noise 
is readily treated in the frequency domain using spectral-density 
functions. Frequency-domain, linear-circuit analysis of equivalent 
circuits greatly facilitates the usual treatments. (Sometimes, as in 
mixer theory, periodic rather than time-invariant situations are con- 
sidered, and frequency-domain analyses are still convenient.) 

In CTD's the situation is quite different. During the transfer of the 
charge from one storage region to the next, the (noiseless) currents, 
charge densities, conductances, etc. associated with the transfer of the 
signal packet are rapidly time-varying. 9 As a result, the noise generated 
during transfer is nonstationary ; that is, the statistics of the current 
fluctuations that give rise to noise vary appreciably with time. Once 
the transfer is complete, a certain amount of noise has been introduced 
into the signal independently of the noise acquired in prior or sub- 
sequent transfers. The processes of interest, therefore, are nonperiodic 
as well as nonstationary. Under such conditions, we have found it best 
to work in the time domain using correlation functions. 

The second basic difference between this treatment and other 
treatments arises from the physical structure and operation of the 
device. Ordinarily, we can externally control the noiseless portion of 
the voltages and currents associated with the signal. In CTD's, 
however, the size of the charge packets, in addition to the charac- 
teristics of the clock voltages, controls the charge transfer. Thus, the 
noise currents that we usually calculate function as additional driving 
terms in calculating charge transfer and, of course, result in fluctuations 
in the sizes of the packets. But it is the reaction of the charge transfer 
to these currents, and not the currents themselves, from which we 
must calculate charge-packet fluctuations. 10,11 (Normally, we control 
either the voltage across, or the current through, the device of interest 
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and calculate fluctuations in the other. Here we control neither 
directly.) 

A third basic difference is that it by no means suffices to calculate 
only the noise introduced during a single transfer. The collecting 7 of 
the noise must also be carefully calculated to include correlation 
effects that lead to a suppression in the cumulative noise added to each 
packet while being transferred by the device. The number of un- 
expected effects makes treating noise in CTD's truly interesting. 12 

At first the nonstationary feature of the noise coupled to the reaction 
of the charge transfer might seem to preclude a reasonably simple, 
analytic treatment that produces useful results. Indeed, some of our 
expressions will be a bit complicated. However, owing to certain 
suppression effects, which will be discussed in detail, the mean-square 
induced fluctuation in the size of a given charge packet at the end 
of a single transfer is nearly independent of the size of the packet. 
Thus, having calculated the magnitude of this fluctuation, it will be 
possible to obtain a simple, meaningful expression for the spectral 
density of the accumulated noise at the output which is independent 
of the signal. This quantity is quite useful for evaluating the effects 
of noise on the analog performance of a CTD. By contrast, for digital- 
performance evaluation, only the accumulated mean-square fluctuation 
in the size of a charge packet at the output is needed. As we shall see, 

this quantity is arrived at in a straightforward manner using our 

time-domain analysis without the necessity of working in the frequency 
domain at all. 

In what follows, we shall outline briefly the noise sources whose 
effect on the output charge packets we shall calculate. Then, following 
a review of the lumped-charge model of a CTD, which has proven to 
be so useful in the analysis and calculation of incomplete transfer 
coefficients, we calculate the statistics of the noise charge introduced 
in a single transfer in terms of the statistics of the microscopic fluctua- 
tions inducing this noise. Thermal, trapping, clock-voltage, emission- 
current, and leakage-current fluctuations are considered. Three 
different types of compounding, storage-process, transfer-process, and 
modulation, are then treated neglecting incomplete-transfer effects. 
Following this, the general problem of compounding in the presence 
of incomplete charge transfer is treated and interesting suppression 
effects are uncovered. The paper ends with calculations of the maxi- 
mum signal-to-noise ratio expected for state-of-the-art devices. As the 
techniques employed in calculating nearly all aspects of the noise are 
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of necessity different from those usually encountered, we have chosen 
to elaborate our methods in some detail. From a knowledge of the 
noise at the output, the operational limitations 13 and error rates 14 
associated with CTD's can be assessed. 

The remainder of this article is long and intricate. This arises not 
from the nature of the noise treated per se, but rather because of the 
complexity of the devices considered. Nonetheless, familiarity with 
the physical operation of CTD's as well as acquaintance with Brownian 
motion and shot noise are the only prerequisites needed. No use is 
made of highly developed, sophisticated techniques of noise theory. 
This is not because such techniques are not considered applicable; 
rather, it is felt that the methods used here are the simplest available 
for a rigorous treatment. In addition, it is felt that these methods will 
be useful in treating noise in other dynamic devices as well. Readers 
well versed in noise theory may feel that this problem is amenable to 
existing methods. While the contrary is not claimed here, we believe 
that such treatments will be more involved than might be expected 
at first sight. 

II. SOURCES OF RANDOM NOISE 

We should at the outset stress that we are concerned here only 
with random noise. We are not concerned with signal distortion arising 
from sampling, incomplete transfer, 15,16 direct clock coupling, nonideal 
regeneration, or any other deterministic process. This is not because 
we feel that such problems are unimportant. Rather, our philosophy 
is that random noise is unavoidable, whereas, in principle, deter- 
ministic "noise" can be greatly reduced or compensated for by careful 
design. 13,14 Thus, it is random noise that plays a major role in setting 
the operational limits of CTD's, which is a problem of general current 
interest. 

Figure 1 indicates schematically the sources of CTD noise with which 
we shall be concerned. At the input, if the charge packet is created 
by photon absorption as in imaging application, or if it is injected by 
an emission-limited mechanism, full shot noise [_((Q — Qo) 2 ) — eQol 
will accompany the signal Q. If, on the other hand, the packet is 
injected into the input via a resistor-controlled circuit, as in regenera- 
tion and storage applications, no shot noise will be introduced. 

The CTD itself simultaneously transfers charge and stores charge. 8 
In its charge-transfer capacity, thermal noise from the Brownian 
motion of the carriers composing the transfer current and trapping 
noise from the fluctuation in occupancy of interface states 17 arise. Care 

1214 THE BELL SYSTEM TECHNICAL JOURNAL, SEPTEMBER 1974 



INPUT 




TRANSFER 










STORAGE 




OUTPUT 


SHOT 




THERMAL 
TRAPPING 




LEAKAGE 
TRAPPING 




CLOCK 
COUPLING 


































COLLECTING 








































MODULATION 












CLOCK-VOLTAGE FLUCTUATIONS 
CAPACITANCE VARIATIONS 









Fig. 1 — Sources of noise for a CTD. 

must be exercised in calculating this transfer noise because of the 
tendency of relaxation effects to suppress thermal and trapping 
fluctuations. 7 Since the transfer of charge from one storage region to 
the next is controlled by the conductance of the storage regions them- 
selves (or portions thereof), shot noise in the transfer current is 
totally suppressed for all practical purposes. (Exception : Should such 
current become barrier-limited, then some shot noise can result. This 
special case is treated in Section 6.1.) In its capacity of charge storage, 
noise from leakage-current fluctuations and trap-occupancy fluctua- 
tions is introduced. Although intrinsically much smaller than transfer- 
noise sources, storage noise is unsuppressed, and, hence, it can be 
important in some cases. At the output, we have the usual problems 
with detector noise, but we shall not consider these problems at 
this time. 

What is of basic interest is, of course, the cumulative noise in the 
charge packet by the time it reaches the output of the device. As it 
turns out, one must be very careful in calculating the collecting of the 
noise introduced into each packet during each phase of each transfer 
cycle. For example, if at the end of a transfer phase, a transfer noise 
of +AQ has been added to the signal, by conservation of charge, a 
quantity of charge — AQ has been added to the charge left behind. 
These two contributions to the noise are highly correlated, and this 
correlation must be taken into account in calculating the noise spectral 
density at the output. 7 Incomplete transfer of charge distorts the 
noise as well as the signal, and must also be included in collecting 
effects. Fluctuations in the clock voltage coupled to fluctuations in the 
sizes of the storage capacitances of the individual storage cells of the 
CTD give rise to modulation noise, which collects quite differently 
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from the noise induced by fluctuations in the clock voltage in the 
absence of cell-parameter variations along the CTD. 

These, then, are the aspects of the noise produced by CTD's that 
we shall discuss. In calculating these various contributions to CTD 
noise, we were surprised by the variety of the results : the presence or 
absence of suppression in individual transfer phases, differences in 
compounding suppression, stationary fluctuations in charge-packet 
fluctuations induced by nonstationary noise sources, etc. Athough we 
by no means treat all aspects of CTD noise, the methods we develop 
should be helpful in calculating the influence of nearly any source 
of noise on the output of a CTD. 

III. MODEL OF CHARGE-TRANSFER DEVICE 

All device noise arises ultimately from fluctuations inherent in the 
transfer of charge carriers between states characterizing the flow of 
charge through the device. 18 It follows that to calculate CTD noise, 
we must first understand how a CTD operates under noiseless con- 
ditions. 10 - 11 As this topic has been the subject matter of a large number 
of papers, 1011 a detailed elaboration here is not necessary. We shall, 
however, briefly review a lumped-charge model 10 ' 11 of charge transfer 
within a single transfer unit of a CTD, which has proven useful in 
discussing and calculating incomplete charge transfer in CTD's. This 
model enables us to express the device current as a function of several 
characteristic voltages, which in turn control the charge passing 
through the device. Using this model, we can then calculate fluctuations 
in the sizes of the transferred charge packets in terms of the current 
fluctuations which accompany the charge transfer and which can be 
calculated by standard means. 18 

The lumped-charge models we shall use to calculate CTD noise are 
shown in Figs. 2a and 2b. Charge Q a stored on capacitor C, is trans- 
ferred to capacitor d during the transfer cycle. (The clock voltage 
which determines the duration of the transfer cycle is V c .) In sub- 
sequent cycles, the charge on Cd is transferred to the right step by 
step to other storage sites, and charge from the left is brought into 
C a — both processes are modeled by repetitions of this model for a 
single transfer unit. For the present, we shall be concerned only with 
a single transfer unit, either of type (a), if the charge transfer is 
characterized by a single-step process, or of type (b), if the transfer 
is characterized by a two-step process. In the latter process, an inter- 
mediate capacitor C& is inserted between C, and Cd to enhance their 
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Fig. 2— Lumped-charge model including the effect of interface states for the 
transfer of charge through a single stage of a CTD. (a) Single-step transfer, (b) 
Two-step transfer. 



mutual isolation so as to reduce the incomplete transfer. As we shall 
see, this also has the effect of reducing the noise acquired by the 
transferred packet. 

Let us now review the dynamic features of the lumped-charge 
model of Fig. 2a in some detail. Stored on C a = C a (V c , V d) V c , V ss ) is 
a time-dependent quantity of stored charge Q, given by 



Q a = [ V 'C B (V s ,V d) V ei V„)dV' s ; 
Jv, 



(la) 



stored on C 8S (V 8 , V d , V C) V as ) is a relatively small, time-dependent 
quantity of interface (trapped) charge Q, a given by 



q,s = C" c aa {v„ v d , v c> v; 3 )dv; s . 

Jv„ 



(lb) 



The total charge Q to be transferred through the conductance 
/ = I(V a , V d) V c , V aa ) is given by 



Q = Q. + Q 8t . 



(lc) 
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Stored on C d is a charge Qd given by 

Qa = Qo- Q = f^ C d (V a , V d , V e , V 8e )dV d , (Id) 

Jv d0 

where Q is the total charge to be transferred through this single 
transfer unit. (Q is some constant for each transfer event: it is the 
quantity of charge initially distributed between Q 8 and Q ss to be 
transferred to C d .) The time dependence of the decay of Q s during 
the charge transfer is governed by the equation 

Qd = -Q = I(V„ V d , V c , V 88 ) + i n (t), (2) 

where i n {t) is the device noise current induced by microscopic fluctua- 
tions within the conductance I. To determine CTD noise we must 
(i) solve eqs. (1) and (2) for fluctuations in Q d induced by i n (t) (and 
by fluctuations in V c and in trap occupancy), and (it) express i n (t) 
in terms of the microscopic fluctuations from which it is induced. 
For two-step transfer processes, as modeled in Fig. 2b, the dynamics 
of the charge transfer are more complicated. One must now be con- 
cerned with Q 8 , Q b , Qd, Q.«i, and Q as2 denned by 

Qs = f V ' Cs(V' s , V b) V c , V 88l )dV' s , (3a) 

J v., 

Q b = f V " C b (V s , V' b , V d , V c , V 8Sl , V 88i )dV' b , (3b) 

Jv b0 

Q d = f Vd C d (V b , V d , V c , V 882 )dV d , (3c) 

JVdo 

Qssi = C' 1 C 88l (V 8 , V b , Vc, V' ssl )dV' ssl , (3d) 

Qss2 = ( Vui C 882 (V b , V d , V c , V' ss2 )dV' ss2 . (3e) 

Jv„ in 



In addition, we have 

Q = (Qs + Qssi) + (Q b + Qss2) m Qx + Q 2 , (3f) 

and 

Qd = Qo - Q = Qo - (Qi + Q 2 ). (3g) 

The corresponding dynamic equations governing the time dependence 
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of Q S) Qb, and Q d are 

Qi = Q a + Q as i - -Ii(7., K», F CJ F„i) - i„i(0, (4a) 

Q 2 = Q b + Q ss2 = Ii(V s , V b , V c , V aal ) + im(t) 

- h(V b , V d) V c , V.. t ) - i n2 (t), (4b) 

Qt = h(V b , V d , V e , V„ 2 ) + i n2 (t), (4c) 

where again the noise currents i„i and i n2 are the device noise currents 
induced by microscopic fluctuations within the conductances I\ and 
h, respectively. To calculate the fluctuations in Q d , we shall solve 
eqs. (3) and (4) for the influence of i n \, i n 2, and fluctuations in V c and 
in trap occupancy Q 8 ,i and Q ss2 , on Q d . 

IV. NOISE INTRODUCED IN ONE CHARGE-TRANSFER CYCLE 

Let us now inquire into the influence of device noise on the size 
Q d of a charge packet transferred through one single-step transfer 
element of a CTD. (Owing to the additional complexity associated 
with a two-step transfer element, we shall relegate the noise treatment 
of this case to Appendix A). It is the purpose of the remainder of this 
section to write the mean-square fluctuation in Q d , q d , in terms of the 
independently fluctuating voltages and currents characterizing the 
transfer and storage process. By expressing all dynamic quantities in 
terms of a minimum number of independent ones at the outset, we 
can greatly reduce the number of cross-correlations which we must 
eventually include. 

If we assume that the fluctuations themselves are sufficiently small, 
then we can linearize eqs. (1) and (2) about their time-dependent, 
noiseless solutions. The form taken by eq. (2) upon expanding to 
lowest order in the fluctuations (and subtracting out the noiseless 
portion) is 

dl .61 . dl . dl . . ,_. 

qd= ~ q = w-, v > + w- d Vd + dv c v ° + w; a v " + *-' (5) 

in which lower-case letters indicate the fluctuation portion of the 
quantity of interest, e.g., Q = Q° + q, Q d = Q d + q d , etc. Since v, and 
v d depend upon q, we must also linearize eq. (1), which will yield this 
dependence. In calculating incomplete charge transfer, it was necessary 
to pay special attention to the dependence of C s and C d on V t , V d , 
V e , V ta . This was because most terms in the small-signal expansion of 
Q t and Q d contributed to the incomplete transfer of charge. For our 
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purposes here, however, we can recognize such terms as small; that 
is, we can drop terms of the order of incomplete transfer times noise 
compared to terms of the order of unity times noise. As incomplete 
transfer per transfer is at most 3 X 10 -3 in devices of interest, this 
is certainly warranted. Proceeding in this manner, we obtain 

q = C.(V°, V° d , V% V°„)v. - C,(V° S0 , V° d , V% V° ss )v so 

+ C(V°, V° d , V°, V° ss )v s , - C, S (V°, V° d , V°, V% )v tto 

= C.v, - C° s v t0 + G at v, s - C°„v tt0 (6a) 

and 

q o = - q = C d (V°, V° d , V°, V%)v d - C d (V% V ^, V° c , V° ss )v do 

= C d v d - C° d v d0 . (6b) 

It is these equations that must be solved for v a and v d in terms of q, 
v c , and v„. 

At this point a minor subtlety enters. In Fig. 2a, contact 1 is tied 
to the clock voltage V c . Since a fluctuation v c in V e cannot instan- 
taneously alter the amount of charge stored on C„ fluctuations v, and 
v„ must be correlated, and indeed this correlation is contained in 
eq. (la). A similar argument also applies to contact 2, v d and v do , and 
eq. (Id). (Since the capacitance C„ is not tied directly to the clock, 
we need not concern ourselves with trapped charge at this point.) 
It follows from eqs. (la) and (Id) that 

°- = l>t v < = W,cl v '= {+l) c]''< (7a) 

and that 

dV do dV d C d , 1X C d / 7 i>. 

"- = W7 v ' = W,cl v ' = (_1) ci v - (7b) 

[In (7b), dV d /dV c =—1 because it is — V c that is connected to 
contact 2.] With respect to the traps, since V tt is the effective level 
to which traps are occupied during the transfer cycles, V„, is un- 
affected by the dynamics of the operation of the device and, hence, 
v..o = 0. We may now solve eq. (6a) for v„, (6b) for v d , and insert 
into eq. (5) to obtain 

-{^-W7.) v " + i " (8a) 
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from which it follows at once that 

q d (t) = |V«p[-|VvT(«")][(«f- + Wl)\*W 

-(^-^l^ t,)+i M' (8b) 

where we have taken g T « g m and denned the time-dependent relaxa- 
tion time t by 

1/ T es g JC B . (8c) 

(The capacitances C, and C d are of comparable size.) In eq. (8), g m 
is the forward conductance (dl/dV.) and g r is the reverse conductance 
(— dI/dV d ). Because of the inherent undirectionality, g m ^>g r . The 
statistical distribution of q d (t), which we seek, can be determined from 
(8b) and the statistical properties of v c , v„, and i n . 

Our expression for the noise fluctuation in the transferred charge 
given in eq. (8) has several interesting features which should be 
carefully noted. The charge fluctuation that accumulates is not simply 
the integral of the noise current over the time interval. Rather, it is 
the integral of the noise current times a damping (relaxation) factor. 
Physically, this arises because if an excess quantity of charge — q( = q d ) 
has been transferred, then the subsequent current (which depends 
most strongly on q s ) is reduced from what it would have been in the 
absence of the fluctuation. The reduced current causes less charge to 
flow from Cs to C d , which in turn tends to partially compensate (null 
out) the effect of the previous fluctuation. This leads to a suppression 
of the device noise, which is shown explicitly in eq. (8b). During the 
initial portion of the charge-transfer cycle, the transfer current is 
largest, as is the conductance g m . Hence, although the larger the 
current, the more noise is present in the charge transfer, because the 
damping is also largest initially, the effect of this noise on the trans- 
ferred charge can be expected to be greatly suppressed. On the other 
hand, near the end of the charge-transfer cycle, while the noise from 
the transfer of charge and its accompanying noise is much reduced in 
size, so is the damping. As it turns out, it is the noise produced during 
the end of the transfer cycle that is most important. This is con- 
venient, because it means that the size of the noise produced by the 
end of the transfer cycle is independent of the size of the initial charge 
packet (as long as some charge is transferred). 

Of primary interest is the mean-square fluctuation in the transferred 
charge (q d (t f Y), where t f is the time at the end of the transfer cycle. 
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Since the fluctuations v C) v, a , and i„ are mutually independent, as we 
shall see when we determine their statistical distribution, it follows that 

(q*(tf) 2 ) = / fl " dh j^ dt 2 exp T - f"d&/r(£))\ 

X exp [ - £ ti/r(ti][( 9m + fi)\ h (g m + §)\ ^ 

x ^<^«>+0-i)L0-4)L 

X (qss(ti)qss(t2)) + <*.(*i)i«(*i)>l (9) 

In (9), we have used C, u v„, = q,, and g m /C, = 1/t. Thus, to determine 
(?d), we must calculate the autocorrelation function of v c , v„, and t». 
For stationary noise, such quantities are well known. 3-6 The purpose 
of the next section will be to calculate these autocorrelation functions 
for the nonstationary conditions that enter the present problem. The 
reader should note carefully at this point, moreover, that the factors 
multiplying these correlation functions involve time-dependent quan- 
tities characteristic of the detailed (but noiseless) solutions to the 
nonlinear device equations [eq. (1)]. Although further simplifications 
can be made in some cases, it should be evident that, in general, one 
must understand the noiseless problem in order to do the problem 
with noise. 

Two other important aspects of our results (8b) for qd(t) and (9) 
for (qd(ts) 2 ) are these. First, as a result of the suppression factor, 
q<i(tf) and, hence, (qd(t/) 2 ) are for all practical purposes independent 
of the size of the signal. Owing to the suppression, qd(t,) depends most 
strongly on details of the charge transfer for t tt tf. However, consider- 
ing the typical size of a, the coefficient of incomplete charge transfer 
(a ^ 10~ 8 ), the details of the charge transfer for t tt t f can deviate by 
only about 10 -8 , which for our purposes is quite insignificant. Second, 
qd(t f ) is a stationary random variable, even though its statistics must 
be derived from the nonstationary distributions of v c , v aa , and t'„. 
These two results greatly simplify our treatment of compounding of 
the noise in Section V. 

(Note that, to avoid undue complication, we have left out leakage 
current into the storage regions, as well as the noise associated with this 
current. As this noise source is uncorrelated with the other sources 
considered above, we can treat it in a separate section. When this 
component of noise is included, it is no longer the case that q d = — q 
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[eq. (6b)], and, hence, greater care must be exercised when we sum 
the noise added with each transfer. This effect also shows up in the 
distinction between transfer process and storage noise, which is 
discussed in more detail in Section V.) 

V. NOISE INTRODUCED BY MICROSCOPIC FLUCTUATIONS 

In the preceding section, we expressed the mean-square fluctuation 
in the charge transferred during a single transfer cycle in terms of the 
autocorrelation functions of the various contributions to the device 
noise. In this section, we shall calculate these correlation functions 
and then use the results to estimate the noise of each type introduced 
during a single transfer cycle. Our task is eased considerably because 
of the extensive effort that has already gone into the study of noise 
in solid-state devices. 3-6 

Although at first sight the character of the charge transfer from 
C 8 to Cd in a CCD appears to be rather different from that in an 
IGFET bucket-brigade device, in fact, the two types of charge transfer 
can be treated in a similar manner. We shall carry over this similarity 
in treating the device noise associated with i n (t): we shall make use 
of the understanding available of noise in IGFET's and then apply 
these results to the CCD as well. In so doing, we must (and shall) 
be careful to make note of certain important differences between 
CCD-mode and BB-mode transfers which can affect the noise 
calculated. 

5.7 Thermal noise 

Of the primary sources of noise present in IGFET's — thermal noise 
at high frequency, 19 generation-recombination (g — r) noise at inter- 
mediate frequency, 20 and 1// noise at low frequency 21 — by far the 
most important source of noise associated with i n (t) is the thermal 
noise arising from the Brownian motion of the charge carriers in the 
inverted region of the semiconductor, which forms the conductance /. 
(This is evident because the spectral densities of 1// and g — r noise 
have dropped considerably from their peak values by 10 5 Hz, the 
lower bound on CTD transfer frequencies.) (In CCD's the proximity 
to the semiconductor-insulator interface of the charge being stored as 
well as transferred greatly enhances the role of interface trapping on 
(g%). In fact, for CCD's, it appears that interface noise is the dominant 
form of device noise in these devices. For reasons that are apparent in 
the next subsection, we shall treat this contribution as a trapping noise 
associated with v aa rather than as 1// noise associated with i n .) 
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Of the early theoretical work on thermal noise in IGFET's, the 
treatment by Jordan and Jordan 22 seemed clearest to us on first 
reading. These authors note that a spontaneous current fluctuation 
i,(x, t) at x along the channel of the IGFET (0 ^ x ^ L) will give 
rise to a voltage fluctuation v a (x, t) = i,(x, t)dx/na(x, t), where a* is 
the mobility of the carriers and <r(x, t) is the mobile charge per unit 
length of channel at time i. The spontaneous voltage fluctuation 
V s (x, t) will in turn induce an z-dependent voltage fluctuation all 
along the channel, which in its turn induces a fluctuation in the 
source-to-drain current id given by 

id(t) =t<t{x, t)v a (x,t). (10a) 

Making use of the above relation between v a and i s , we obtain 

id(t) = i a (x, t)dx/L. (10b) 

[The same result can be obtained using the impedance field method 
(IFM). 23 Alternatively, one can develop a current-current method 
analogous to the current-voltage method employed in the IFM. This 
is outlined briefly in Appendix B.] The contribution to the noise 
current i n (t) due to id(t) induced by fluctuations all along the channel 
is from (10b) 

i n (t) = f is(x,t)dx/L, (11a) 

Jo 
so that 

(in(ti)i n (ti)) = I dxi I dx2(i a (xi,ti)is(x 2 ,t2))/L 2 . (lib) 
Jo Jo 

The autocorrelation function of i s (x, t) for thermal noise can be 
found from that of the current density j s (x, t) obtained from micro- 
scopic noise theory :8 and given by 

0'.(xx, h)j t (x2, U)) = 2kTn P (x h ti)8(x 1 - x 2 )5(^ - h), (12) 

where p(xi) is the charge density at Xi. Inserting (12) into (lib) 
and noting that p(x) = <r(x)/A(x), A(x) being the cross-sectional area 
of the channel at x, we obtain 

(in(h)i n (t»)) = 2kT8(h - h) f »<r(x, h)dx/L\ (13a) 

Jo 

The integral in (13) must be found from a knowledge of the noiseless 
operation of the device at time h. In general, the integral is directly 
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proportional to the forward conductance 22 so that 

(in(h)i n {U)) = 2kT5(h - U)g m (h)H n {h), (13b) 

where H n = 1 for the IGFET operated in the linear region, H n = § 
for the IGFET operated in saturation, and H n = \ for diffusion- 
limited current. 11 In most cases of interest, the CTD is operated in 
saturation so that E n = f. For very long clock periods, H n = £, 
appropriate near the end of the cycle, can be used. As the 30 percent 
uncertainty between \ and f is tolerable at present for noise calcu- 
lations in these devices, we shall not concern ourselves with the 
additional dependence of H n on the time t\. (Setting H n = § and 
H n = \ gives adequate upper and lower bounds on the thermal noise.) 

[Although (13b) was derived for the bucket-brigade type of charge 
transfer, we shall assume that it is valid for CCD-type transfer as 
well, where, of course, one uses g m appropriate to the CCD device. 10,11 
This assumption clearly breaks down, however, when (<$)* becomes 
comparable to or larger than the amount of free charge incompletely 
transferred in a single transfer phase. The reason for this difficulty 
goes back to the assumption of linearity in eq. (5). This, of course, 
also applies, mutatis mutandis, to all other noise contributions as well. 
Treating these interesting nonlinear problems is, unfortunately, 
beyond the scope of this paper.] 

Returning now to (9), we can calculate the contribution to (g}) 
arising from thermal noise during a single transfer cycle. If we set 
Hn = h ignore the time dependence of C, (a 0.1 percent effect if 
a = 10 -8 ), and recall that t _1 = g m /C„ inserting (13b) into (9) and 
integrating over ti and t 2 , we obtain 

{&)\kT = IkTC. (14a) 

in the limit that 



expT- f" dt'/rit')] «1, 



as it must if the incomplete transfer a is to satisfy a <3C 1. (In deriving 
the above, we have not assumed that g m is independent of time.) A 
corresponding result was first given by Boonstra and Sangster; 24 
however, their ($) was four times larger than the right-hand side of 
(14a). We believe that (14a) is in fact correct. 

[If one ignores the nonstationarity of the noise, a result similar to 
that of Boonstra and Sangster can be obtained in the following manner. 
For thermal noise, the mean-square current fluctuation is given by 
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(id) = 4:kT%g m B, where g m is the conductance and B is the bandwidth. 
If r = C,/g m is the characteristic relaxation time, then B = t _1 and 
(<$) — (»1)t 2 = 4/cTf C s . (The quantities (gjj) and (id) are mean-square 
fluctuations, not Fourier coefficients. The quantity (i\) is obtained 
from the spectral density of id using the Wiener-Khinchin theorem.) 
While this result is not too bad, attempting to calculate noise in more 
complicated, nonstationary situations using stationary results can 
lead to trouble. For example, for a two-step transfer process, we obtain 
for the thermal contribution to (g}) the result [from eq. (62)] 

M)\w = lhTC h + i*rc.(l + 6- 1 )- 1 , (14b) 

where b = {g m 2/Cb)(g m i/C,)~ l , and is assumed to be independent of 
time. This result is more of a challenge to obtain from arguments 
based on stationary noise sources.] 

5.2 Interface (trapping) noise 

We now focus attention on calculating the statistics of 

q„(t) = C„(t)v„(t), 

the fluctuations in occupancy of the interface states during charge 
transfer. To proceed, we must first write down the dynamic equations 
governing the trapping. These are simple since we are working in the 
lumped-charge approximation. If greater accuracy is desired, one can 
be more microscopic and include the position and energy dependence 
of the trapping states, as well as their capture cross section and 
thermal-release time. The procedure is similar but much more involved. 
The dynamic equation relating the flow of charge into and out of 
the traps is 

Q„ = It(I, Q„) - I r (Q„) + in,,, (15) 

where the trapping current I t is given by 

I t = I(<r x /A)(Q\ S - Q..)/<&, (16a) 

and where the release current I r is given by 

I, = Q„/r T (t). (16b) 

In (15), i n „ is the noise current associated with the filling and emptying 
of traps; in (16a), / is the transfer current, <r x /A is the ratio of the 
effective cross section of all the traps to the area A through which I 
flows, and the factor (Q l ss — Q„)/Qf n is the fraction of the total number 
of active sites Q' ss that are empty; in (16b), r r is the thermal-release 
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time. Writing Q aa = Q° ss + g u , expanding to first order in the noise 
q s „ and subtracting off the zero-order equation, we obtain for the noise 

dI t dl . dl t dl r . m v 

-[-7& + 4 + ^]*" + l - (17b) 

— «H/r«(0+U ^ 17c ) 

where r, is the effective trapping time (dI t /dQ„), and where riJ 1 
equals the bracketed quantity in (17b). (Ordinarily I T « /, so that 
the first term in r,! 1 can be neglected. Under these circumstances, the 
noise contribution from different types of traps will be independent, 
and the traps may be treated independently.) 

Proceeding now, we solve (17c) for q S s(t), obtaining 

q„(t) = J* dt' exp \-j' tf"/T..(f")l ins.it'). (18) 

To obtain the correlation function of q„(t), we must calculate 
(q„(ti)q„(tt))- The resulting expression involves the correlation func- 
tion of i nB ,(t), which, since i nS8 is an elementary microscopic process, 
is given by 18 

<<«.(*iX...(fa)> = eZWt) + J?(«0]8(«i " ti)i d9) 

where e is the charge on an elementary carrier. (The above is not 
obtained from a spectral density, since for a time-dependent current, 
a nonstationary process, spectral density is not defined.) For t 2 > h, 
it follows that 

(q..(h)q..(tt)) = exp [- j* dt/r t ,(t)\ 

■ j ll dt exp I" - j^ dt'/r sg (t')] eU°t(t) + l?(0J (20) 

If this result is inserted into (9), summation over all trapping states 
yields the contribution to the transfer noise arising from these states. 
As is evident from (20), in general, one must know the noiseless 
solution, especially t„(0, W) and I° r (t), in order to actually calculate 
the noise. A few simplifications, however, permit us to recover 
Tompsett's result 17 for a single-step, CCD-mode transfer, which is 
valid in the limit of I t « / n i.e., near the end of a transfer cycle. 
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To calculate the contribution to (<$) arising from interface-state 
trapping, we must, according to (9), know not only (Qii(h)q„(t%)), 
which we found above (20), but also (dI/dQ 8a ). This quantity entered 
into a previous discussion 10 of the role of surface states on incomplete 
transfer. There it was noted that one could write 

dl/dQ.. =-a(t)/ T (t), (21) 

where for CCD-mode transfer a = and for bucket-brigade-mode 
transfer a falls from a = C s /C C h during most of the transfer cycle 
(C C h being the channel capacity) to a £^ toward the end of the cycle, 
where the IGFET current becomes emission limited. [The vanishing 
of a for CCD-mode operation is due to the fact that the field felt by 
a mobile carrier due to another carrier is independent of whether the 
other carrier is free or trapped. The large value of a for BB-mode 
transfer arises from the fact that relatively small changes in the 
threshold voltage, induced by changes in trap occupancy in the channel 
region, can result in relatively large changes in transfer current. If a 
is large, it enters (9) as a 2 , while, as we shall see, (g? s ) is proportional 
to C. («C e A for the IGFET channel). This implies that the ratio (R) 
of the surface-state-noise contribution to (g§) for bucket-brigade to 
that for CCD is C,/C C h- If, on the other hand, the BB-mode is not 
turned off until the channel current becomes emission limited, then 
a ^ 0. If we assume that the suppression factors in (9) damp out the 
noise introduced while a ^> 1, then the above ratio (R) becomes 
Cch/C„ which means that the contribution to (q 2 ) of q aa can be ignored 
for BB-mode transfer. We shall assume that this is the primary 
operational region of interest, pointing out, however, that if the 
BB-transfer mode is terminated while a is large, one should expect an 
increase in that portion of ($) arising from trapping.] 

Returning to CCD-mode transfer where we expect interface states 
to play the largest role, if we assume in eq. (9) that (q aa (ti)q lt (t2)) 
varies in time slowly compared to r(t), then we can integrate over 
t\ and U to obtain 

(qMf) 2 ) W (q,.(t f y). (22) 

If, in addition, we focus attention on long transfer cycles (10~ 6 second 
or longer), taking 1° » I", and r t » r r , then from (20) we obtain 

(q..(tf) 2 ) =e f"dt 

•exp [-2 J*' <k'/T r (t')] Q° ss (t)/Tr(t) + (q.,(t c y). (23) 
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If we make one more assumption that at time t c (0 < t c < tf), 
Q»» = Q°n = Qu t and for subsequent t(t c < t < t f ), Q° ss decays as if 
It = 0, then 

QUt) = QS. exp [ - f* dt'/r T {t') ] , (24a) 

(qM 2 ) = 0. (24b) 

Inserting these into (23), we find that 

(Qssitf) 2 ) 

= qQ l s S exp\- j'' dt'/Tr(t')M 1 - exp [ - J" <fc'/r r (o]}, (24c) 

as found by Tompsett 17 using other reasoning. Should conditions be 
such that any of these assumptions are unwarranted, then, of course, 
the result [(20) in (9)] is more complicated. 

The result given in (24c) has the following significance. If on the 
average r r » (t f - t c ) or r r « (t f - t c ), then ($,) for such states is 
negligible. In the former case, the transfer process occurs too rapidly 
for the traps to respond; in the latter case, the trap occupancy can 
follow the transfer current quite closely; in either case, the noise is 
greatly suppressed as a consequence. Thus, only when r r ~ tj — t c 
can the traps influence the charge transfer. 

If we now sum over the distribution of traps, assumed uniform in 
energy E, and take r r _1 proportional to exp ( — E/kT), then following 
Strain 25 we obtain 

<<&) = e*kTN t ,A In 2, (24d) 

Tompsett's result, 17 where iV™ is the number of interface states per 
unit energy per unit area, and A is the active trapping area. Like 
Tompsett, 17 we conclude that the noise introduced into Q d from 
trapping has a mean-square value of 

(ql) = e*kTN,.A In 2, (25) 

that is, 

5.3 Clock-voltage noise 

The influence of fluctuations in clock voltage on the noise added to 
the transferred charge in a single transfer can be treated very quickly, 
especially if it is white as we shall assume. If the spectral density of 
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the noise at zero frequency is S v (0), then 18 

(v c (ti)v c (t 2 )) = |5„(0)5(ii - t 2 ). (26) 

From (9) it follows that the contribution to (<$) from clock- voltage 
fluctuations is 

<g2)|oiook = *£.(0) [*di 

Jo 

•exp [-2 J" dt'/T(t')] (g m + dI/dV e )*\ t . (27a) 

We recall that g m = CJr. The direct dependence of / on V depends 
strongly on the type of device, and one must be careful not to include 
in dI/dV c terms already included in (dI/dV,)(dVJdV c ). Thus, for a 
single-step bucket brigade, dI/dV c =+g m , while for a single-step 
CCD, dI/dV c = 0. If we set dI/dV c = bg m , then 



dt 

o 



<<z§>|ciock = *s.(0)C!(i + by f l 

•exp [" -2 j" dt'/T{t')\ T-*(t). (27b) 

In general, (27b) must be evaluated from a knowledge of r(t) obtained 
from the noiseless charge-transfer characteristic. If, however, we 
assume that we can write dI/dQ a = dI/dQ„ then 

dt/r(t) = g m dt/C s = dtdI/dV a /C a = dtdl/dQ. = dtdl/dQ, =-dI/I. 

It follows that 

<^>|oiock = £»S,(0)C 2 ,(1 + by f'' j^jj^fy (27c) 

where // = I(t = t f ) and I = I(t = 0). In this form it is clear how 
the integral yields the effective bandwidth B of the white clock- voltage 
noise. One can replace the t~ 1 in (27c) by dI/dQ 8 , which in turn may 
be calculated as a function of i" using eq. (8) of Ref. 11 (Q, = C a V 8 ). 
For our purposes, we shall be content with 

(Si)loiook = tf,(0)BC!a + by. (27d) 

Of course, if the clock-voltage noise is not white, (26) should be 
replaced with the actual correlation function, which then introduces 
an additional time-dependent function in the integrand of (27). 
Relation (27d) will be useful in discussing modulation noise. 
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VI. OTHER NOISE SOURCES 

In the preceding section, we discussed the three primary sources of 
noise that contribute to the fluctuations in the size of a single charge 
packet as the result of a single transfer. In this section we examine 
two other sources of noise which do not readily fit into (9) without 
undue complication. For simplicity, we shall not include the thermal, 
trapping, or clock-voltage fluctuations discussed in Section V, as these 
can be superimposed linearly on the results discussed below. 

6.7 Emission-limited-current shot noise 

In describing the types of noise expected to be generated in CTD's, 
we noted that shot noise would not play a significant role under 
ordinary circumstances. This is because, in simplest terms, a single 
transfer unit of a CTD is like two capacitors connected by a resistor, 
and, owing to strict charge neutrality in the resistor, shot noise is not 
present. If, however, the clock period U of an IGFET CTD becomes 
long (t > 10 -6 second), the channel current becomes partially emission 
limited at the source towards the end of the transfer cycle. 11 In this 
case the shot noise associated with the emitted current will not be 
totally suppressed. We shall now show how this can be treated. At 
the same time it will become clear why "ideal" resistors totally sup- 
press shot noise. 

In Fig. 3 we represent the barrier region between the diffused source 
and the channel by a conductance with current I e and the channel 
by a conductance with current /. The voltage V, at the source end of 
the channel is set V a as in Ref. 11. In Section V we ignored I e , assuming 
its conductance was much greater than that of /. We now consider 
the more realistic situation in which the conductances of I c and / are 
comparable near the end of the transfer cycle. (In such circumstances, 
the g m introduced in Section IV and used extensively in Section V, must 
be replaced by the series conductance of I e and / in the expressions 
for thermal, trapping, and clock- voltage noise.) 



I 






Ic<V p .V a .V c ) KV ai V d ,V c ) 

Fig. 3 — Lumped-charge model for a CTD including a barrier conductance /«. 
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What is happening physically can be outlined as follows. Voltages 
V a as well as V p and Vd can fluctuate in response to fluctuations 
induced in the conductances. If a small, spontaneous current fluctua- 
tion i, (positive, say) occurs through I e , then V a will increase, simul- 
taneously inducing a (positive) fluctuation i in i" and a (negative) 
fluctuation i e in I e . The net current fluctuation is i = i, + ie- If the 
conductance of I e is much larger than that of /, then V a will change 
very little, ifSi 0, and hence i e fa—i,. In other words, a current i t is 
induced by i a which nearly cancels i, ; with the larger conductance of 
I e , a small v a induces an i e sufficient to suppress the i, fluctuation, 
while with the smaller conductance of /, v a induces a much smaller i 
and /. In the limit of an "ideal" resistor, since it is the bulk which 
controls the current and not the contacts, i — * in the limit of zero 
contact resistance, and the shot noise is totally suppressed. 

Let us now calculate the shot noise introduced during emission- 
limited operation. We start as before with the equations of the model : 



Ie(V p , Va) = Ie(V p - V a ), 


(28a) 


I = I(V a , V d) V e ), 


(28b) 


Q. = [ V *c.(v p ,v ai v c ), 

Jv po 


(28c) 


-& - I e + i. = /• 


(28d) 



In (28a) we have assumed that I e is primarily a function of (V p — V a ), 
anticipating that for barrier current, I e is well-described by a diode 
equation. Since the variations of / with Vd or V c and of C. with V a 
or V c are small, we shall ignore these dependences as we did in Section 
IV. Setting I t = I° e + i e and / - 1° + i, it follows from (28) that 

dl 

U = gy- (v P ~ v a ), (29a) 

i = ^- v a , (29b) 

q. = C,»„ (29c) 

-q. = ie + i. = i. (29d) 

Solving for q, in terms of i„ we obtain 

-a = *L (dI/dV P )(q./C p )+i s 
q> dV a dI/dV a + dI/dV p ' K<sua,) 
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or letting g m = dI/dV a , g b = dI/dV p and solving for q„ 

q* = -«. = / o ' dt' exp \- fa *7n(f")l C*.(l + w/^^lf, (30b) 

where r 6 = [({7m 1 + Qb l )C p ~]. This expression tells us the fluctuation 
induced in the transferred charge packet qa resulting from the shot 
noise associated with the current emitted into the channel from the 
source. In the gr& — > oc limit (zero barrier resistance), qa — * 0. In the 
gb — > limit (infinite barrier resistance), <?<*—» f i,dt, full, unsuppressed 
shot noise. 

In deriving (30b) our primary goal was to determine ($)■ Since 
(i.(ti)i,(h)) = el°(ti)8(ti - h), 

(gi)lshot = / o " dt exp r -2 J*' dt'/M?)] el°(t) 

■ll + g b (t)/g m (t)l- 2 . (31) 

To evaluate (31) requires, in general, a knowledge of the noiseless 
solution to the charge-transfer equations. We can obtain a feel for the 
result if we recall van der Ziel's expression for the current spectral 
density of diode shot noise : 6 

Si(f) = 2e(I° + 7.) = 2kTg b J^rff - (32a) 

(I is the diode leakage current that we have neglected.) Although 
derived for stationary noise, it is straightforward to redo the derivation 
for nonstationary noise to obtain an autocorrelation function of 
the form 

(t.(*i)».(**)> = e[I-(0 + /.)S(ti - h] 

- kTg b (h) T^q^f *(«i - *») (32b) 

or 

el°(t)8(h - t 2 ) « kTg b (ti)8(tx - h), (32c) 

where in (32c) we have (again) ignored the diode leakage current 
7 . Now then, given the barrier, (g§) will be largest if g b « g m . Inserting 
(32c) into (31) and integrating on t, it follows that 

(g})\ B ^<ikTC p , (33) 

which is comparable with the thermal noise produced by I (%hTC p ). 
If one may assume that towards the end of the cycle g b (t) = dg m (t), 
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then 

<^>Uot = hkTC p (l + d)~K (34) 

While (33) and (34) should be sufficient for estimating the contribution 
to Qd of shot noise, if more accuracy is desired (31) should be used. 

6.2 Leakage-current noise 

Although leakage currents are small and, hence, fluctuations in 
them even smaller, it is of interest to briefly analyze leakage-current 
noise. This is because (i) a portion of the contribution to g<j from 
leakage-current fluctuations is not suppressed and (ii) no longer does 
qd = — q s . We shall restrict ourselves to the role of leakage current 
which flows into the source, drain, and channel of an IGFET, bucket- 
brigade, single-step transfer CTD. Leakage current into the source 
and sink of a CCD can probably be treated in a similar manner. 

If we assume that a leakage current per unit length of the channel 
J C h(x, t) enters the channel, following Jordan and Jordan, 22 we find 
that a fraction (L — x)/L flows towards the source and x/L flows 
towards the drain. We have assumed, of course, that f J C h{x, t)dx <K /. 
From the standpoint of noise this means that 

q a = -q./r + / dxj ch {x, t)(L — x)/L + / dxj,{x, t), (35a) 



and 



q* = ^+f dx 2%^* + f dxjd(x, t), (35b) 

T Jo Li Jd 

where J = J° + j. Solving (35a) for q, and (35b) for qd, we obtain 
q.(t) = £ dt' exp T - /J dt"A(*")l 

• f J q L dx ^-^ j ch (x, t') + J[ dxj x (x, Ol (35c) 



and 



q*(t) = /J dt' [^j + j Q L dx ^^ + j d dxj(x, o] • (35d) 



Comparing (35c) and (35d), we note that while the fluctuation j is 
suppressed in q, and, therefore, in the first contribution to qd, it is 
not suppressed in the second or third contributions to qd. Thus, even 
though leakage currents may in themselves be small, since a portion 
of their noise is not suppressed, one may expect to see some contribution 
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to qa from this source. As a rough estimate, if after many transfers 
the size of Qd is increased by Q e (leakage charge), one may expect 
(ql) due to leakage to be about qQ e . The autocorrelation and cross- 
correlation functions of qd and q a can be calculated at once from (35b, c) 
if one notes that for full shot noise, 

(j(*ii h)j(x 2 , U)) = eJ(x h t!)8(ti — i 2 )5(xi — x 2 ). (36) 

Since an accurate calculation of the results again requires some 
detailed knowledge of the noiseless charge transfer and the role of 
leakage current, we shall not pursue this topic further. 

VII. STORAGE-PROCESS AND TRANSFER-PROCESS NOISE 

Up to this point in our discussion, we have been concerned solely 
with the various contributions to qd that result from a single charge 
transfer of a single charge packet. A charge packet reaching the 
output of a CTD has, however, been transferred typically 10 2 to 10 3 
times. The noise in the output packet is an accumulation not only of 
the noise acquired by the packet of interest during each transfer, but 
also the noise contained in incompletely transferred portions of pre- 
ceding packets. In addition, there is correlation between the noise in 
successive packets at the output. Some of this correlation arises, of 
course, from the incompletely transferred portions picked up along the 
line. 7 However, even in the absence of incomplete charge transfer, 
there is substantial correlation from packet to packet (see Fig. 4). For 
example, for thermal and interface-state noise, we noted that through- 
out a single transfer cycle qd = — q- While qd accompanies the packet 
of interest, q{ = — qd) is picked up by the next packet. 8 At the output the 
correlation will affect the spectral density of the total noise accompany- 
ing the signal. In this section, we consider this correlation and in the next 
section, we consider modulation noise, both in the absence of incom- 
plete transfer; in Section IX, we discuss the output noise including 
incomplete transfer effects. (Finally we remind the reader that all 
along we have been concerned with random noise which is generated 
in addition to the signal distortion resulting from incomplete transfer. 
In many cases, physical processes which contribute to random noise 
also contribute to incomplete transfer. However, it should be kept in 
mind that while the noise is random and characterized by stochastic 
processes, the incomplete transfer is deterministic and characterized 
by a specific transfer function for the entire device.) 

We noted in the case of contributions to qd arising from thermal, 
trapping, clock-voltage and barrier-current fluctuations that conserva- 
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Fig. 4 — Schematic illustrating the difference between the origin of storage-process 
and transfer-process noise. 

tion of charge led to the relation (fa + g = 0, where q* is the noise 
added to the transferred charge packet and q, which is picked up by 
the following packet, is the noise added to the untransferred charge. 
Using the notation of Ref. 8 in which qff is the transfer-noise charge 
introduced during the o-th phase of the sth transfer cycle as charge 
flows from the (/* — l)th to the /xth phase of the mth CTD cell, the 
accumulated noise AQ? in the charge at the last phase of the iVth cell 
at the end of the rth transfer cycle is given by 

AQT=E E (gW*— m ~ q?-\ N - m) -iJ. (37) 

m=l fi=l 

[In (37), p is the number of charge transfers per cell. If n is the total 
number of charge transfers in the CTD, then n = Np. During the 
(r + l)th transfer cycle, Q?( = Q? + &Q?) flows from the device to 
the output through a low-pass filter. If by g(t) we designate the output 
current per total charge transferred for a single packet, then the 
output noise current is given by 

i(t) = E AQfy(< - rto), (38) 

r 

where t„ is the clock period. One can calculate the spectral density of 
i(t), obtaining 

Si(f)\r P = 4m,(q*)Tvf \gtf)\*(l ~ cos 2tt/// ), (39a) 

where TP denotes transfer process and f is the clock frequency 
(fo — to~ 1 )- [Here we have assumed that (q 2 ) is independent of 
(s, 0-, m, n) or, equivalently, of the signal. This is quite reasonable 
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since the suppression effect renders the noise dependent only on the 
final portions of the charge transfer. And owing to the smallness of 
the coefficient of incomplete transfer, the dependence of the final 
portion of charge transfer on the signal is negligible for noise purposes.] 
For / <JC f /2, Sir? is greatly suppressed below its mean value, >S,(/ /4) 
(see Fig. 5). This effect results from the fact that at low frequencies 
one is averaging over such long times that nearly all the q<i are cancelled 
by their corresponding q =—q<i. For /£d f„/2 (f is constrained by 
the relation | /| < / /2), jS»tp is in fact enhanced by the strong, 
mutual correlation between adjacent packets. The suppression of the 
transfer process noise at low frequencies is advantageous, since by 
increasing f one can increase the signal-to-noise ratio (S/N). 

For the moment, we ignore the transfer-process noise associated 
with leakage and consider only storage-process (SP) noise, which we 
define as fluctuations that influence the size of each packet inde- 
pendently during each transfer phase. Since, under such circumstances, 
there is no correlation between the noise in different packets, the 
spectral density of the filtered current at the output is independent 
of frequency (white) and is given by 



S,(/)|sP = 2n(< ? 2 > S p/o|0(/) 



(39b) 



l\g(f)\ « 1 for | /| < /o/2]. If <g 2 ) T p = (<Z 2 )sp, the integral of &(/) 
over < / < fo/2 is twice as large for TP as for SP noise. This is 
because in the case of TP noise each fluctuation contributes to two 
charge packets, whereas in the SP case, only one packet is affected. 
We note that £,(/) for SP noise is not suppressed for /« f„/2. Thus, 
although leakage noise is expected to be small, since neither is a portion 
of it suppressed in forming (g 2 )sp nor is £,(/) suppressed by packet- 
packet correlation, the role of leakage-current noise may in some cases 
be more important than is usually appreciated. 
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Fig. 5 — Noise spectral densities Sbp and <Stp plotted versus frequency /. 
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In discussing TP and SP noise above (and in Ref. 8), we have 
assumed that different, or rather independent, mechanisms contribute 
to (q 2 )sp and (<? 2 )tp, and, for TP noise, that qd ■= — q. For leakage- 
current noise, we found that in fact the same fluctuation could con- 
tribute to both SP- and TP-like noise, and that qd 7^—q = — q,. The 
question is raised : What is the spectral density of such noise? 

To calculate the spectral density of the noise in cases where TP- 
and SP-like effects are correlated, we proceed as follows. Let the noise 
accumulated on the source in a single transfer phase be q,(s, a, to, n) 
and the noise on the drain be qd(s, a, to, h), where (s, a, m, fi) is denned 
above. Then 

N p 

AQr V = L. L {qd[r - (n - to), /z, to, mJ 

+ ff.[r- (N-m)- 1, m, m, /*]}, (40) 

where we continue to ignore incomplete transfer effects. It follows 
that the spectral density of i(t) denned by (38) is given by 

SiU) = 2nf \g(f)mq 2 d ) + (q*) + 2(q d q a ) cos (2tt/// )]. (41a) 

FJFor pure TP noise, q a = — qj, and we recover (39a) ; for pure SP 
noise, q„ « qd owing to suppression, and we recover (39b).] By rewrit- 
ing (41a) slightly, it separates into SP and TP portions : 

8t(f) = 2nf \gU)\*{(qd + q*Y) 

+ 4n/.| fl i(/)|»<-g.g i )[l - cos (2r///.)]. (41b) 

Thus, we find that even when the sources of noise leading to TP and 
SP noise are correlated and, hence, more complicated than those 
treated in Ref. 8, we still obtain an unsuppressed, white-noise con- 
tribution and a suppressed contribution. If incomplete transfer effects 
are included, the spectral densities become much more complicated. 
We take up this matter again in Section IX. 

VIII. MODULATION NOISE 

In Section 5.3, we calculated the noise introduced into the signal 
during a single transfer resulting from fluctuations in the clock voltage. 
To calculate the effect that this modulation noise has on the output of 
the CTD, we must notice that each clock-voltage fluctuation is felt 
simultaneously by each transferring packet throughout the entire 
device. If the elements of the CTD are physically quite similar, each 
fluctuation will induce nearly the same noise contribution to the 
(flVi Qd — —Qa) pair generated in each CTD element during each cycle. 
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This means that modulation noise can be expected to be substantially 
suppressed, which is in fact what we find. In calculating modulation 
noise, we shall first assume that the transfer parameters of each pair 
of transfer elements of the CTD are identical, and then we shall 
include some variation in these parameters. 

If by q we represent that portion of qa introduced by clock-voltage 
fluctuations, since q is purely TP noise, it follows from (37) that 

AQ? - t \ -flfc2*r* + q?t + E* (#%-«>., ~ flP±Jf-«.jl (42a) 

H~l I n=l J 

= i fagr - q^nJ. (42b) 

M-l 

The simplification attained in (42b) results from the above observation 
concerning the similarity of contributions to the noise throughout the 
CTD, which stated quantitatively is that 

q?f = gtt (43) 

for all cells ?n, n during a given cycle s. It follows from (42b) that the 
contribution to the total mean-square fluctuation in the output due 
to clock- voltage fluctuations is simply 

<(AQ?) 2 >| clook = 2p(^)| c i ock , (44) 

where (ql) is given by (27d) and p is the number of phases per cycle. 
The most interesting feature of (44) is that it is independent of n, 
the number of transfers each charge packet is subjected to in the CTD. 
Thus, this portion of modulation noise is not compounded and as 
such can be expected to be small. 

If we introduce the possibility that the physical parameters charac- 
terizing each transfer stage of the CTD are slightly different, then 
(43) will not hold, and, as a consequence, the clock-voltage fluctuations 
will not essentially cancel. Let us assume that these physical param- 
eters are distributed according to some probability distribution. Then, 
in place of q in (42a), let us write (q + q') where q is the charge fluctua- 
tion averaged over the distribution of the device parameters, and q' 
is the deviation from this average. Inserting (q -f- q') into (42a) for 
q, expressing AQ as (AQ + AQ'), a straightforward calculation leads to 

(AQP6Q?) = 2p($ 2 )[5 r . s - !(«,.._* + «r..+Ar)] (45a) 

and 

(AQ; n AQ' s n ) = 2n(q~ r2 )l8 r ,s - K*r..-i + «r,.+i)l (45b) 
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in which the bar d enotes a veraging over the device parameters. To 
calculate (cp 2 ) = ((q - q) 2 ) = (q 2 ) - (q 2 ) for clock-voltage fluctua- 
tions, it is perhaps simplest to carry out the average first over the 
clock- voltage fluctuations obtaining (27), and then average over the 
device parameters. If we assume that the primary variation from cell 
to cell is the size of the storage capacitance C„ then (27d) yields 

{{AQ' r N ) 2 ) = 2n-hS v (0)B(l + 6) 2 <(C 8 - C,) 2 ). (45c) 

This result expresses the effect on the size of the charge packet of the 
coupling between fluctuations in clock voltage and deviations in 
storage capacitance. Unlike ((AQ?) 2 ) given in (44), {{AQ' T N ) 2 ) is 
proportional to n, expressing the fact that during each transfer event, 
noise is added to the signal. Finally, we note from (45b) that AQ' T S has 
the character of transfer process noises so that its spectral density is 
suppressed for / « / /2. The spectral density of AQ?, which can be 
obtained from (45a), is weird. It oscillates once between and 4p(§ 2 ) 
each time / changes by only f /N. It is probably best approximated 
as white of size 2p(q 2 ). 

IX. INFLUENCE OF INCOMPLETE CHARGE TRANSFER ON COLLECTING 

In the two preceding sections we have ignored the influence of 
incomplete charge transfer on the noise accumulated in charge packets 
by the time they reach the output. In Section VII, we took a = 
to illustrate as simply as possible the suppression in the spectral 
density of transfer-process noise at low frequencies. In Section VIII, 
we took a = to treat the multicorrelated charge fluctuations induced 
by the clock voltage in a straightforward manner. We now include 
incomplete transfer and find that, even though it alters only a small 
fraction of the signal and, hence, of the noise on each transfer, typically 
10~ 3 to 10~ 4 , its accumulated effects are extremely important in 
some cases. 7 

The effects of incomplete transfer on a charge packet are not simple, 
even when linearity is assumed. The key is to write down a general 
expression for the charge, seen as a function of time at the output, 
which arises from a given charge packet introduced at an earlier time 
on an arbitrary cell. Once this is done, one can combine the effects on 
the output of fluctuation-induced noise charge which is created on 
each stage during each phase of each cycle. The result is a complex, 
combinatorial expression which for na ^ 1 is nontrivial to evaluate 
numerically. Despite these complications, the results obtained are 
worth the effort needed to obtain them. 7 
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To calculate the effect of incomplete charge transfer on the collecting 
of noise, we first calculate the portion of the noise q? added to the 
charge packet in the rath cell during the rth clock cycle which is 
observed at the output during the sth clock cycle. As defined, q? is 
related to the q?f of Section VII according to 

q? = f 8^. (46) 

Thus, for purposes of simplicity, we shall ignore for the moment the 
fact that the transfer of charge within each cell involves p distinct 
(independent) transfers. The error involved in so doing is of the order 
of na 2 , which for devices of interest is less than 0.01. Since the linear 
model we shall use to treat incomplete transfer is probably not this 
accurate, this approximation is justified. However, we must be careful 
to use (46) when we calculate terms such as (q?q n T ) so that we do not 
neglect correlation effects. 

The dynamic equation that governs the transfer of the charge 
q(r, to) on the rath cell during the rth cycle is 

q(r + 1, to + 1) = q(r, to) — eq(r, to) + eq(r, to + 1), (47) 

where e is the coefficient of incomplete transfer per cell. 15 [One can 
relate e to a, the coefficient of incomplete transfer per transfer, 16 
through e = pa, or more accurately, through (1 — e) = (1 — a) p -] 
Using (47), it is straightforward to calculate the charge q m (s, N + 1) 
observed during the sth cycle at the output, the (N -f l)th stage of 
an iV-cell register, as a function of charges q (r, m), which are added 
to the packet present in the mth. cell during the rth cycle. The result is 

q m (8,N+l) = (1 - e ) Ar ~ m+1 r 5 ( , ' + ^~ W ) 

• e'qds - (r + N - to + 1), to]. (48) 

[If e = 0, the limit of negligible incomplete transfer q m (s, N + 1) 
becomes 

q m (s, N + 1) = g [s - (N - m + 1), to J 

Thus, the additional output charge seen during the sth cycle is just 
the charge added to the with cell, N — m + 1 cycles earlier.] Using 
(48), we calculate the total noise charge AQf observed at the output 
of an iV-cell CTD during the sth cycle by replacing q (r, m) by the 
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noise charge q? and summing overall cells m. The result is 
AQT = (1 - •) £ (1 - •)*- ± ( T + N r ~ M ) 

m=\ r=0 \ T J 

■e'q?- {r+N _ m+lh (49) 

where 

q? = qd(r, m) + q.(r - 1, m). (50) 

In (50), qd(r, m) is the noise added to the charge on drain during the 
cycle of interest, and q s (r — 1, m) is the noise added to the charge 
remaining on the source during the preceding cj r cle and picked up by 
the charge packet coming by on the next cycle. [The expression in 
(49), as well as our treatment of collecting, is valid whether or not the 
statistics of the individual q/s and q s 's must be calculated using a 
nonlinear approach.] From eq. (49) we can calculate nearly all com- 
pounding effects of interest, a few of which we now consider in some 
detail. 

For digital purposes, the most important quantity of interest is 
<(AQf ) 2 ), the mean-square fluctuation in the size of the output charge 
packet. This quantity can be calculated from (49) keeping in mind 
that the only nonzero, cross-correlation function that enters is 
(qd(r, m)q a (r, m)) for all r and m. The result is 

<(AQ") 2 > = <(<M + qs) 2 )H s Ap, N) + 2(-q s q d )H TP (p, N), (51a) 

where qa and q s are noise added to the charge in the drain and source, 
respectively, during a single transfer cycle, the statistics of which we 
calculated in Section V and Section VI, #sp is the collecting factor 
for storage-process noise q a = 0, and Hj P is that for transfer-process 
noise ( — q, = qa). The analytical expressions for H SP and #tp, 
ignoring an unimportant factor of (1 — e) 2 , are 

n "° / r 4- N — m\* 

HsAP, N) =p E (1 - e) 2( "- m) £ ( J * 2r (51b) 

m=l r=0\ T / 

and 

#tp(p, N) - V £ (1 - e) 2 <"—> 

m = l 

■i.( r+ r m H i - f+1 r + + r m )- ™ 

These can be evaluated exactly (Appendix C) as well as approximately 
(Appendix D). [The latter is necessary because the former, although 
exact, is difficult to evaluate for the large N and small e of greatest 
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interest. These calculations are nontrivial. For example, except when 
Ne « 1, attempting to calculate (51b) or (51c) using Stirling's formula 
to approximate the factorials is doomed to failure.] One may note at 
the outset that for e = 0, Hs p = #tp = pN = n the number of 
transfers, as expected. This result is also approximately valid for 
Ne = na « 1. 

The approximate results obtained from evaluating the sums in 
(51b, c), 

H 3P (p, N) = P x N (e) = p(N + h)t- l e~ b Uo(b) + h(b)2 - \ (52a) 

and 



H T p(p, N) = P<pn(c) 






= p{(l - e) + r l (x + l)-^[X Ar _ 1 (e) + 1 - 


- Xw(e)]}, 


(52b) 


where 






t m (1 - «) 

(i + *y 




(52c) 


b = 2(N + |) € (1 - e)- 2 , 




(52d) 


and 






*~£ + :l r=* 




(52e) 



are quite interesting. (X N and ips are given in Appendix D.) For 
Ne <<C 1, where we expect incomplete transfer to play a very minor 
role in the compounding of noise, indeed we find that H S p and Hip 
are nearly equal to pN = n, the total number of transfers experienced 
by each packet in the device. This is just what one expects: The 
cumulative, mean-square noise charge after n independent transfers is 
just n times the mean-square noise charge following a single transfer. 
[The factor of two in the TP term of (51a) arises because, as explained 
in Section VII, for each +q noise contribution there is a — q contribu- 
tion. Thus for each transfer, two noise terms are produced.] As N 
increases (e is fixed), however, incomplete transfer plays a more 
significant role, altering the noise in two important ways. First, of 
course, the noise is incompletely transferred along with the signal. 
For Ne 2> 1, H S p increases only as n*, reflecting this attenuation of 
the noise. Second, and even more important, for TP noise, incomplete 
transfer enables each -\-q, — q pair created during each transfer of each 
packet to mix and, hence, null out or suppress the total noise. Thus, 
for Ne ^> 1, H T p approaches constant value, (2a) -1 , independent of N. 
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In other words, for Ne = na X> 1, the collecting saturates and further 
collecting is totally suppressed (see Fig 6). 7 Of course, the signal is 
greatly distorted by the incomplete transfer if Ne ^>> 1. Nonetheless, a 
CTD has maximum storage capacity 13 for 1 < Ne < 4, and by using 
dynamic detection or optimum linear filtering (transversal filtering), 
one can greatly suppress signal distortion from incomplete transfer. 13 
Thus, calculating the compounding factors, H S p and H T p, for Ne 
other than Ne <3C 1 is not an academic exercise. 

For analog purposes it is necessary to calculate the autocorrelation 
function (AQf AQ; V ) using (49) and (50), from which the current (or 
voltage) spectral density of the noise can be obtained as in Section 
VII. We shall pursue this no further than to point out that since 
q„ = — q must be transferred one more time than q, the effect of in- 
complete transfer on each q, q s pair will be slightly different, and this 
will reduce their mutual correlation at the output. Thus, total sup- 
pression at zero frequency is no longer expected. 

Incomplete charge transfer will also affect fluctuations at the 
output caused by modulation noise. In Section VIII we found that 
modulation noise was so highly correlated that, in the absence of 
incomplete transfer, the largest portion of modulation noise was not 
compounded. Introducing incomplete transfer, however, will destroy 
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Fig. 6 — Suppressed, collecting factors H S p(n) and #Tp(n) plotted versus n, the 
number of transfers, for a = 10 -3 . 
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the exact correlation which led to this cancellation. Again, if desired, 
the effect of the collecting of incomplete charge transfer on the modula- 
tion noise for digital or analog applications can be calculated using 
(49) and (50). 

Noise present in the input signal can clearly be treated simply as 
part of the signal with respect to incomplete-transfer effects. A some- 
what different situation arises, however, when the signal enters the 
CTD all along its extent, as is the case in imaging applications. Here 
full (photon) shot noise results in a signal-dependent contribution to 
the mean-square, input noise of each packet equal to eQl m , where 
Q° s _ m is the average size of the packet injected at site m during cycle s. 
However, in this case, the input packets undergo different numbers of 
transfers, p(N — m), between the cell in which they are injected and 
the CTD output. Thus, the influence of incomplete transfer affects 
the signal and noise originating at each different cell differently. As 
with the other compounding effects, the collecting of this shot noise 
in the presence of incomplete transfer can be calculated if desired 
using (49) and (50). 



X. CALCULATIONS 

In the foregoing we have discussed a number of sources of noise in 
CTD's and their effect on the output signal in the absence and presence 
of incomplete charge transfer. In most cases, inserting the appropriate 
physical parameters for the noiseless operation of the device of interest 
suffices to calculate (q 2 ), the mean-square fluctuation in size of a 
given charge packet acquired on a given transfer. Then, using (51), 
the influence of this fluctuation on the output signal can be determined. 
Such calculations are, in general, difficult, owing to the necessity of 
evaluating integrals such as those in (27b). However, realistic approxi- 
mations can be made as indicated to obtain useful results. 

It is of interest, however, to determine the minimum amount of 
noise expected to be present in CTD's assuming one can minimize 
clock-voltage fluctuations, surface states, incomplete transfer, etc., 
and operate each device so as to avoid emission-limited currents, etc. 
In this ideal situation, one is left only with thermal noise, or with 
thermal noise plus shot noise on the input signal, and incomplete 
transfer (intrinsic and modulation in the sense of Refs. 10 and 11). 
The signal-to-noise ratio (S/N) was calculated for two characteristic 
coefficients of incomplete transfer ( = 10~ 3 , 10~ 4 ) and four characteristic 
capacitances (C = 1, 0.1, 0.01, 0.001 pF) as a function of the number 
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of transfers n from input to output : 

(i) Using eq. (14a) for the thermal noise acquired per transfer. 

(it) Using eq. (51a) for the influence of this noise on the output. 

(Hi) Taking the thermal noise to be purely TP, including shot 

noise at the input when present. 
(iv) Including incomplete transfer effects on both the signal and 
the noise. 

The results are plotted in Figs. 7 and 8. The ratio Q /C designates the 
maximum signal level (10 volts), and one-half this amount (5 volts) 
is the minimum signal level. Since the square of the signal charge is 
proportional to C 2 , while the mean-square of the noise charge is 
proportional to C for both thermal and input shot noise, the S/N 
decreases proportionately with smaller C (small CTD cells). As the 
number of charge transfers is increased, the contribution of device 
noise to the total noise soon dominates that of the input noise. (From 
S/N one can also calculate the maximum information storage capacity 
of the CTD as a function of n, a, and C.) In general, other noise 
sources are present which reduce S/N from the ideal results shown 
here. While valid in general for BB-mode transfers, in the case of 
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Fig. 7 — Signal-to-noise ratio for CTD (a = 10 -3 ) with storage capacitance of 
1 pF, 0.1 pF, 0.01 pF, 0.001 pF. 



1246 THE BELL SYSTEM TECHNICAL JOURNAL, SEPTEMBER 1974 




10 3 10" 

NUMBER OF TRANSFERS 



1CP 



Fig. 8 — Signal-to-noise ratio for CTD (a = 10 -4 ) with storage capacitance of 1 
pF, 0.1 pF, 0.01 pF, 0.001 pF. 

CCD-mode transfers, if the rms thermal noise becomes comparable 
to the free charge incompletely transferred on a single transfer, then 
the noise predicted by our linear approach will overestimate the 
true noise. 

XI. CONCLUSIONS 

In the foregoing we have calculated the influence of several im- 
portant sources of noise on the output of a CTD in a manner which 
includes all important relaxation, suppression, and correlation effects. 
In so doing, we encountered interesting effects, such as nonstationary 
noise, and we uncovered a number of unexpected results, such as the 
nearly total suppression of the spectral density of transfer-process 
noise at low frequencies and the total suppression of the collecting of 
transfer-process noise after many transfers of a charge packet. Because 
of the novelty of these and other effects, they were treated in con- 
siderable detail in a manner which did not presuppose considerable 
prior experience with formal treatments of device noise. In noting the 
results, the proportionate decrease of the S/N with decreasing storage 
capacitance was illustrated. This unavoidable feature will ultimately 
limit the practical size of CTD's. 
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Several of our results, despite the relatively complicated manner in 
which they were obtained, are remarkably simple, e.g., thermal noise, 
(q 2 ) cc kTC, surface-state noise, (q 2 ) « ekTN as A, etc. One wonders 
whether there exists a general approach which circumvents the 
necessity of paying such careful attention to the detailed processes 
that accompany the transfer of charge through a CTD. We feel that 
there does not exist such an approach. Nonetheless, in Appendix E 
we outline a very rough method which ignores nearly all details of the 
charge transfer event. We find the intuitively appealing result that 
the mean-square noise acquired per transfer is to within a factor of 
the order of unity equal to full shot noise on the (total and not differ- 
ential) incompletely transferred portion of the charge. Although 
providing a rough rule of thumb, since the approach is not totally 
reliable for calculating incomplete transfer, its accuracy for treating 
noise is not guaranteed. By contrast, the methods used in the bulk 
of this paper should be applicable in many types of integrated-circuit, 
dynamic devices of which CTD's are the first examples. 

There are several interesting noise problems that we did not con- 
sider here. For example, in discussing modulation noise we indicated 
the possibility of nonuniformity in the physical parameters of each 
cell coupling to clock-voltage fluctuations to produce a collecting 
source of noise. Such cell nonuniformities also, of course, will result 
in a distribution of a's, the coefficients of incomplete charge transfer. 
This will, in turn, result in an additional effect on the nature of the 
compounding of the noise acquired in each transfer. The results are 
expected to be no less surprising than the effect such a distribution of 
a's has on the signal. A distribution of a's about their mean a leads 
to less signal distortion than if all the a's were a . 16 ' 26 (This is actually 
not too significant for application purposes, since usually all deviations 
from the desired a = a<j will be to larger a, thereby increasing a and 
enhancing distortion.) 

A second problem worthy of attention is how to treat the noise 
in cases where the clock-voltage waveform does not turn off the flow 
of transferring charge abruptly. 9 In such cases, the nonlinear terms in 
the noise fluctuations are not small relative to the linear terms. The 
noise problem is then nonlinear and much more complicated. 

A third problem, straightforward but tedious, is to calculate the 
S/N's, error rates, and device storage capacities for devices including 
regenerators, optimum linear filtering, and/or dynamic detection. 13 ' 14 
Considering the many new features that have arisen in the present 
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study, as well as the many types of noise that enter in such different 
ways, we feel such problems will not prove unrewarding. 
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APPENDIX A 

Noise Introduced During a Single, Two-Step, Transfer Cycle 

The purpose of this appendix is to calculate the noise fluctuations 
induced in the transferred charge during a single, two-step, transfer 
cycle. The procedure followed is the same as that in Section IV of 
the text for a single-step process. Here we shall merely outline the 
essential steps leading to the desired result. 

We begin by linearizing the dynamic equations (4a, b, c) to obtain 

/ dli , dl x J dh . dh \ 



(53a) 
qz =-id- qi, (53b) 

Wl Vi + W-/' + Wl"' + W^ v "') + l " 1 - (53c) 

The constraint equations (3a-e), when linearized and when terms of 
the order of the incomplete transfer are dropped, yield the following 
relations among the fluctuating quantities : 

q. = CsV 3 - C° a v, , (54a) 

q b = C b v b - C° b v b0 , (54b) 

q d = C d v d - C° d v d0 , (54c) 

q ss i = C„iU M i - C°ssiV° ssl , (54d) 

q S s2 — CaaiVasi — C%2Vss2- (54e) 

In addition, by the same reasoning that led to eqs. (7a, b) in the text, 
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we find that 

dV so BV s C a ,, n C 8 ,„ . 

v so = -Tfy v c = ^y- g- Vc = (+1) qo Vc, (55a) 

dV bo , dV b C b , C b , ..... 

Vb0 = w; v < = -dv c ci Vc = a ci v <> (55b) 

dV d0 dV d C d , n c„ ... . 

v d0 = -^y- v c = -^y- -q v c = (-1) ^5 v e , (55c) 

where a = dV b /dV c . [For a two-step process, transfer cells whose 
second step is characterized by a CCD-type of charge transfer, a = — 1 
and v'c = v c . For such cells whose second step is BB-like, the value of 
a and v' c will depend on the details of the device. For example, for the 
C4D, a = ; f or the tetrode BB, a = 1 and v' c = y 6 * is the voltage on 
the tetrode bias line. Thus, for purposes of generality we must use 
the form given in (55b). J Also v° ssi and v° s2 can be set equal to zero, 
as discussed in the text. 

The next step is to solve eqs. (54a-e) for v 8 , v b , etc. and insert the 
results into (53a-c). The result of this substitution is 

0ml / . dli\ 

qi= -c7 qi -\ gml + dv c ) v ° 

+ f g m i -gr ~ d y l ) y *«i — *'ni> (56a) 
92 =-q d - qi, (56b) 



<7m2 , ' - dlz ., 

qd = c; q2 + gm2aVc + *r/« 



+ (^ 2 CT " W^) Vs82 + *" 2 ' (56c) 



in which we have assumed that the forward conductances 

flmi = dh/dV s and g m2 = dI 2 /dV b 

greatly exceed the reverse conductances 

g r \ = — 8Ii/dV b and g T i=— dI 2 /dV d . 

It is now straightforward to solve (56a) for qi(t), insert into (56b) to 
obtain q 2 (t), and then solve (56c) for q d (t). The result is 

q d (t) = fc dt' exp [ - j t [ d*"/r.(t")][ ^ (- qi ) + (g m2 a + ^ ) v e 



I / ' \ / C S s2 dl 2 \ i-l 

+ gm2a{V c — V c ) — lg m2 -g -Ty— J V ss2 + *n2 



1250 THE BELL SYSTEM TECHNICAL JOURNAL, SEPTEMBER 1974 



(57a) 



where 



-gi(0 = f* dt" exp 1-^'' dr7[n(r')]j[(g B i + w) Vc 

~( ?ml % 1 "afe) y " 1 + inl ]|^' (57b) 

l/ri = g m i/C, and 1/t 2 = & m2 /C 6 . (57c) 

If we define the two suppression factors according to 

Si(M') = exp(~Jj*7ri(n) 



and 



let 



S 2 (M') = exp(- j t t dt"/r 2 (t")), 

n(t, V) = J] dt" exp (- |' [l/r a (*'") - l/r»(* w '):]*'")/n(l w ), 
insert (57b) into (57a), and regroup terms, we obtain 
q d (t) = ^ (tt ; Is^t, *')»(«, OMO + Sift t')i»»(t') 

-S 2 (M')(^- 2 -^-)| »...(0 

\ Ti aV SS 2/\f 

+ \ [ S,{t,t')n{t,t'){g ml + ^\ ti 



(58a) 
(58b) 

(58c) 



(59) 



Although it is most convenient to keep the Sin product as two terms, 
their combination is of interest : 

Si(t, t')n(t, t') = [' dt" exp f - C dt"'/T 1 (t"') 

- J^dt'"/T2(t"')\/r2(t"). (60) 
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It is clear that to calculate {<$), it is necessary to know (ini(ti)i»i(U)), 
(in2(ti)ins(t2)), (v M i(h)v„i(t s )) t (v SS 2(ti)v SS 2(t 2 ) >, (v c (k)v c (t 2 )) , (v' c (ti)v' c (t 2 )) , 
and (v' c (ti)v c (t2)). These are discussed in the text. 

It is of interest to calculate the contribution to (<$) that results 
from thermal noise. From (59) and (13b) (H n = §), we obtain 

<d)l*T = 2kT | J" dt[SUt f , t)n 2 (tf, t)g ml (t) + 8 2 »(t fi t)g m2 (tn (61) 
Now ordinarily a two-step device is constructed so that 

g m \/C, = TT 1 < T 2 _1 = Qml/Cb. 

Let us assume that the decay current is such that rf x (0& = Ta l (f), 
where b can be taken independent of time (b > 1). Then we can 
perform the integrals in (61), obtaining 

<fl5>| *r = ihTC b + ffcTC.(l + 6- 1 )- 1 (62) 

for the thermal contribution to the noise acquired during a single, 
two-step transfer. The first term is expected from our result (14a) 
for an individual, single-step transfer. However, the second term 
includes a suppression factor whose presence one certainly would not 
have expected a priori using arguments assuming stationarity. 

APPENDIX B 

Outline and Modification of the Impedance-Field Method 

As originally presented, 23 the impedance-field method (IFM) con- 
sists of dividing the problem of calculating device noise into two 
simpler problems. First, it was recognized that a given fluctuation 
in the velocity of a charge carrier at a given location induces a calcu- 
lable fluctuation in the open-circuit voltage. Second, it was shown 
how, from a knowledge of the spectral density of the velocity fluctua- 
tions of the individual charge carriers, the spectral density of the 
open-circuit voltage fluctuations could be calculated. The first problem 
requires only an understanding of the operation of the device of 
interest under noiseless conditions, while the second can be deduced 
from the microscopic behavior of charge carriers in a small region of 
the device. 18 Thus, as long as the microscopic noise is simple, which it 
nearly always is, the device-noise problem is reduced to integrating 
the microscopic noise sources over the device, weighted by the in- 
fluence of a unit fluctuation in each volume element on the output 
noise. 
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Let us now outline the IFM in more detail. 23 If u t is the time-de- 
pendent fluctuation of carrier velocity v< from its mean v *, 

u« = v t - v ot , 

then the dipole current 8p a in a small volume element AV a produced 
by the carriers in this volume is given by 

5? a = q L U,«, (63) 

i 

where the ./-sum is over all carriers in AV a . To calculate the effect 
of 8p a on the open-circuit voltage of the device, we note that if 81 a 
is injected into the device at i a (and if the device is grounded at 
r = 0), the voltage induced at the contact labeled AT, 8V N , is given by 

8V N = Z Na 8l a . (64) 

If Z Na is evaluated for all a, one maps out the "impedance field." If 
now one injects 81 a at r + 5r and removes 81 a at r a , linear super- 
position of small signals implies that 

8V N = lZ Na (i a + 5r) - Z Na (T a )~]8l a (65a) 

= (vZ Nt ).8r8l a = vZ Nt .8K (65b) 

The last equality follows because the dipole current 8 P a equals 8l a 8r 
if 81 a is chosen appropriately. Since we can relate 5f>„ to 8Vn, from 
a knowledge of the statistics of the former, we can calculate those of 
the latter. In particular, from (63) and (65b) we can immediately 
write down the spectral density Sv N (f) of 8V N in terms of that of u„: 

Sv N (f) = T. (V i Z A r )(V i Z A r a )*g 2 n a ,S aUl . tty (/), (66) 

et.i.j 

where i, j each run over x, y, z, and n a is the number of carriers in 
AV a . Since in the text we work in the time domain, we work directly 
with (65). 

In calculating CTD noise, it is most convenient to work with 
short-circuit current fluctuations rather than with open-circuit voltage 
fluctuations. One may redo the above, calculating the short-circuit 
current fluctuation 5/jv induced by 81 a . The result is that 

8I N = B Na 8l a , (67a) 

where 

B Na = Z Na /Z N , (67b) 

and where Zn is the impedance of the device between contact N and 
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ground. In terms of dipole currents, one obtains 

81* = vB Nl -8r8l a = vB Nl -8? a . (68) 

For an IGFET, 

vB Nt = L-H, (69) 

where L is the channel length and / is a unit vector along the channel. 22 
In complete analogy to (66) one finds that 

Si N (f) = E (ViBi, a )(VjBs.)*tfn a 8 mt u i (f). (70) 

In arriving at (66) and (70), we have taken into account the inde- 
pendence between spontaneous fluctuations which occur in separate 
regions AV a . Thus, the expression for (uityjujfa) ) for the IFM, 
which corresponds to (12) in the text, is that for thermal noise 

{Ui{t x )Ui{h)) = =■ M«(*l ~ h)8 {j , (71) 

which lacks the spatial delta function. The equivalence between the 
Langevin method used in the text and the IFM outlined there is dis- 
cussed in some detail in Refs. 27 and 28. 

APPENDIX C 

Exact Calculation of Hsi and W T p* 

In this appendix, we shall evaluate the following sums exactly: 

As the first step, set n = N — ?n, y = e 2 , and perform the sum on r. 
This yields 

*" (e) - n ? ^rr- 37- ( (i - y) -» ) (73a) 

., 1 - 6 , V (1 - e) 2 " d" / 7 " \ 

wW -r^ + £ -in— d7- i (i - 7)- +i ; 

-<$((T^> C»> 

* Derivation due to N. S. Thornber. 
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The second step involves taking the n derivatives with respect to 7 
using the relation 29 



d n , . y t ( n \ d 'f dn ~ l 
dy» U9) to\s) dtfdy*- 



and recognizing that the resulting sum on s is to within certain 
simple factors equal to a Legendre polynomial, 29 P n (x), where 
* - (1 + 7)/(l - r). This yields 

x„(e) = (I-7)- 1 2?*»P.(«) (74a) 

n=0 

w(e) = (1 - 7)- 1 [(1 - «) + EVP.(«) - e«"^=^ P:(x)l , (74b) 

where « = (1 - e)/(l + e). In (74b), if we note that 30 
K(x)/n = [xP n {x) - P^i(*)]/(rf - 1), 

then all that remains is to evaluate Yint n P n (x). 

Before proceeding, we should call attention to a potential source 
of trouble. Since < 7 = e 2 < 1, it follows that x > 1 and t < 1. 
While Legendre polynomials for \x\ > 1 are well-defined, their 
properties are not nearly so simple as they are for x in the usual region 
of interest, |o;| ^ 1. Thus, for fixed e (and 7), if n becomes large, 
evaluating P n (x) for x only slightly larger than 1 is quite tricky. This 
provides motivation for the approach adopted in Appendix D. 

Returning to the remaining sum in (74), we take the generating 
function 31 for Legendre polynomials, valid for \i\ < 1 and \y\ < 1, 

£ t'P.(y) = (1 -2yt + t*)-\ 
s=o 

multiply both sides by P m (y), and integrate on y from —1 to +1. 
Using the orthogonality of these polynomials, we obtain 

«■ = ^^ f\ dyPm(y)(l - 2yt + i 2 )-*. (75) 

Using the summation formula 32 for Legendre's polynomials, it follows 
that 

W-l r\ 7VT/9 

E t n Pn(x) = / dy{\ - 2yt + « 2 )-* -^- 
n=o J-i x — y 

. [Pjr(s)Pjr-i(y) - P N (y)PN-i(x)l. (76) 
To perform the final integral over y we make use of (75), noting that 
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2tx = 1 + t 2 , and hence that 

(1 - 2yt + t*)*(x - y) = 2»(x - y)W. 
Thus, since 

£ dy(i - y)-»P.(y) = -2 ^ £ dy(x - y)-*P n (v) 

and since from (75) 

2t*»(2m + l)- 1 = f 1 <&/[>(* - y)J-*P«(y), 

it follows after some algebra that 

ZjW.fc) = ^72 [^(x)F - P K -i(x)t^J (77) 

71=0 1 — I 

Using (77) for the sums present in (55a, b) yields exact expressions 
for X N and <pn in terms of two or three Legendre polynomials, respec- 
tively. The difficulties encountered in evaluating these expressions 
for N 2> 1 and x > 1 made it clear that another form of the result 
was needed, one in which N and e enter on an equal footing, preferably 
as a product. Such a result is derived in Appendix D. 

APPENDIX D 

Approximate Calculation of H SP and H T p 

In this appendix, we evaluate the sum 

*Ar-i(e) = ?<•?,(*), (78) 

n=0 

where, as in Appendix C, t = (1 — e)/(l + e)> & = (1 + t)/(1 — y), 
and 7 = « 2 . The form we obtain will be a good approximation for N » 1 
and will be very easy to evaluate numerically. 
If we define \p-i = 0, then we can write 

+ n - *«_i = t n Pn{x) (79) 

and evaluate the z-transform rp z of ^ n defined by 

tz = £ ^„2" n . (80) 

n=0 

Thus, from (79) we obtain 



>.-i*.-£„(J)™=[>-H+(i)7 
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(81) 



Solving for \f/ 2 , we obtain, using 2xt = 1 + t 2 , 

+ z = (1 - r-»)-»(l - ftsr 1 )-*. (82) 

Thus, 

}pn = ^if +* zn ~ ldz > (83a) 

where the closed contour in the complex plane includes 2 = 0. Letting 
z = e°, (83a) becomes (a = + ) 

*n - 7T-- [* ** dse«'(l - e-)-*(l - ftr*) - * (83b) 

an exact expression for \lw(e) if n = iV. 

An exact integration of (83b) would recover (77). However, a very 
useful approximation for large n can be obtained at once if we note 
that under such circumstances most of the contribution to the integral 
comes in the vicinity of s = 0. Thus, expanding exp ( — s) « 1 — s and 
then, taking the limits of integration to be from (<r — *<* ) to (<r + i • ), 
yields 



+ n 



^-. f "* ds ««Hr»r« + (r 2 - 1)"H. (84a) 

Z7rt J a-ix, 



This inverse Laplace transform can be found in the tables. 33 The 
result is 

rK « ^ne-^/^ft) + Ji(6)l (84b) 

where 

6 = n (t~ 2 - l)/2 = 2ne(l - e)" 2 . (84c) 

We can match this result, valid for n » 1, to the rc « 1 limit to obtain 
an approximate result good for all n. At this point, however, it must 
be stressed that the major difficulty in approximating \f/„ for large n 
arises from the necessity to cancel exactly any exponential dependence 
of \p n on n. Clearly, physically compounding can increase with n no 
more rapidly than linearly. This delicate cancellation is in fact achieved 
in (84b) since I (b) and J x (b) go as exp (&)/(2tt&)-* for b » 1. There 
are other ways of arriving at (84), but this is one of the simplest. 

To determine the form of \f/ n for n sufficiently small (ne « 1), that 
is, such that an expansion of ^„ in powers of ne rapidly converges in 
a few terms we note from (71), (73a), and (78) that we have the relation 

*jr-i = (1 - 7) I* i\d- e)"e'( r t U )T (85a) 

n=0 r=0 . \ 1 J \ 

= Nil - t(N - 1) + 0(iV 2 e 2 )]. (85b) 
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On the other hand, expanding (84b) to similar order yields 

*„ « n{\ + 2e - en). (86) 

Since (84b) and consequently (86) can be expected to be valid only 
for ne » 1, we are not required to set n = N — 1 in (84b). [Nonethe- 
less, it is encouraging to note that even for ne « 1, if we do 
set n = N — 1 in (86), then (85b) and (86) agree as to the coefficients 
of N and N 2 e, the two largest terms.] Even further improvement can 
be obtained if we replace \f/ n by \f/ n + 8i and n by N — 1 + 82, where 
61 and 5 2 are of the order of unity. This alteration clearly does not 
alter significantly our result valid for large n (?ic» 1), and it improves 
the result for ne <$C 1. After a little algebra, one obtains 81 = — \ and 
52 = f, which ensure that terms of order unity and Ne match. We 
conclude that x N (e) is given quite well by 

x N (e) = (N + h)t- l e-"lh(b) + /i(6)] - |, (87a) 

where 

b = 2(N + i)e(l - e) 2 and t = (1 - e)/(l + c), (87b) 

in which we have ignored the (1 — 7) -1 factor in (73a). (Recall that 
7 = e 2 tt 10~ 6 .) The expression (87a) enables us to calculate the 
compounding factor for storage-process noise for arbitrary N and e 
with relative ease. 

The compounding factor for transfer-process noise is obtained from 
<pw(e), which can be expressed in terms of X N (e) in the form 

<p N (e) = (1 - e) + lX N (e) - 1] 

- ey~ l (x + l)-MaC*Jr(«) ~ 1] ~ ««r-i(«)] (88a) 

= (1 - e) + e~ l (x + l)-H\:X N ^(e) + 1 - X„(e)], (88b) 

where we have used (73) and ignored the pref actor of (1 — 7) -1 . Since 
7 = t 2 , the prefactor of the last term of (88) is approximately (2e) -1 . 
Nonetheless, this term is well-behaved even where our approximate 
result (87) is used in (88) to obtain an approximate <pN(e). 

It is of interest to calculate X^ and <pN in the limit of iVe » 1. The 
former limit is just AT in both cases. In the latter limit 

X N {e) -> (N/re)* = p- l (n/ira)i, Ne » 1 (89a) 

and 

<PN(e) -> (2*)- 1 = (2pa)" 1 J Ne » 1. (89b) 

The second limit is most easily obtained by noting that x(x + l) -1 
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= (1 + 7 )/2, (x + l)- 1 = (1 - 7 )/2, «y-»[(l + 7) - «d - 7)3/2 = 1, 
and (x N _, + 1 - Xat) = 1 + O(AM), We » 1. 

APPENDIX E 

Another Approach to CTD Noise 

Throughout our treatment we have paid careful attention to the 
details of the transfer and storage of charge in CTD's in order to 
calculate the noise induced by the device. We have noted, however, 
that by making reasonable approximations, we can often reduce 
rather complicated, exact results to much simpler, approximate 
results adequate for most purposes. One cannot help but wonder, 
therefore, if there does not exist an approximate but adequate short 
cut for calculating CTD noise. One clue is to compare Tompsett's 
results 17 for incompletely transferred charge Q, with his results for 
mean-square noise charge induced by interface states (q 2 ). If we 
ignore factors of the order unity we find that 

<<Z 2 > = eQi, (90) 

where 

Qi = ekTN aa . (91) 

In other words, the noise induced by interface states is just the shot 
noise associated with the incompletely transferred portion of the 
interface charge. 34 While one will never recover thermal noise from 
such arguments, one expects not to do too badly in estimating the 
influence on the output of those noise sources, such as interface states 
or emission-limited (barrier-limited) currents, which are closely tied 
to incomplete charge transfer. Anyway, the basic idea, that of shot 
noise on incompletely transferred charge, is appealing, however ap- 
proximate and/or incomplete it may in fact actually be. 

We can derive such a result as follows. We make the assumption 
that we can approximate the current / which flows from one storage 
region to the next during a single transfer as I(Q), where Q is the 
charge to be transferred at any given time during the transfer phase. 
(If this assumption had been made in calculating the coefficient of 
incomplete charge transfer, 10 - 11 then certain of the results obtained 
would have been erroneous. Nonetheless, without this assumption, 
details of the charge transfer enter, which we wish to avoid.) As in the 
text we also assume that we can linearize the equations governing the 
noise. Thus, writing Q = Q° + q, the equation of motion 

Q=-I(Q) (92) 
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becomes 

Q0 = _ /( Q0) = _/o ( 93a ) 

for the noiseless transfer, and for the noise 

P-|«+ d q (t), (93b) 

where d q (t) is the naturally arising statistical driving term. Since / 
depends only on Q, the derivative appearing in (93b) is total and is 
evaluated using Q°(t) from (93a). Solving for q(t) is straightforward: 

q(t) = [' dt'd q {t') exp ( - j\ (dI/dQ°) t »dt"\ (94a) 

in which we have used (93a) in going from (94a) to (94b). The 
[1 (!) / 1 (t')~\ factor suppresses the shot-like noise associated with the 
transfer. 

The statistics of q follow from those of d q using eq. (94b). The 
statistics of d q are such that 

{dSMM) = e/[Q°(<i):*(<i - fa). (95) 

(This is not an additional assumption, but rather (95) follows from 
our initial assumption that / depends only on Q.) It follows that at 
the end of the transfer cycle 

/•QO 
o°dQU (Q o i)/P(Q)J, (96) 

where Qo is the initial charge Q°(0) to be transferred, and Qi is the 
mean charge Q°(t f ) left behind at the end of the transfer. [We note 
that in the absence of the expression factor, (q 2 ) = e(Q — Qi) ~ eQo, 
full shot noise. 3 

To proceed, we must know I(Q°). If, toward the end of the transfer 
cycle, Q — > Qi ^> Qd, where Qd is the charge packet size above which 
the primary force driving the transfer current arises from the packet 
itself, then I(Q) will be proportional to Q 2 , and, using (97), 
(q 2 ) ^ eQi/3. If, on the other hand, toward the end of the transfer 
cycle, Q — > Qi « Qd, so that the primary force driving the transfer 
current is diffusion or fringing fields, the I(Q) will be proportional to 
Q, and, using (96), (q 2 ) tt eQi. Taking a more specific example, let 
us set CV = Q in eq. (8) of Ref. 11, as a realistic approximation. It 
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follows that 

<<Z 2 > - e I dQ—) 2/^V ' (97) 

Equation (97) can be integrated exactly. For our purposes, it 
suffices to consider two limiting cases : Qi » 2kTC/e, in which case 
(q 2 ) « eQi/3, and Qi « 2kTC/e, in which case (q 2 ) « eQi. Thus, we 
find in fact that (g 2 ) can be viewed roughly as the shot noise on the 
incompletely transferred charge. However, it should be noted that Qi 
is the total, and not the much smaller differential (aQ ), charge in- 
completely transferred. While the results of this appendix are appealing 
as a short cut, the reader is strongly advised to keep the basic assump- 
tion [7 = I(Q)2 firmly in mind and to use extreme caution in generaliz- 
ing this approach to other problems. 
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